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PREFACE 


The  analysis  described  in  this  report  was  conducted  at  the  US  Army 
Engineer  Waterways  Experiment  Station  (WES).  The  data  described  herein 
were  collected  during  March  and  April  1978  by  a  joint  effort  of  the 
US  Army  Engineer  District,  Portland  (NPP) ,  and  WES.  The  US  Coast  Guard 
provided  a  vessel  and  crew  for  an  installation  of  velocity  meters.  The 
work  was  funded  by  NPP  as  part  of  the  ongoing  WES  model  studies  of  the 
Columbia  River. 

Personnel  of  the  Hydraulics  Laboratory  of  WES  performed  this  study 
during  the  period  1979  through  February  1980  under  the  direction  of 
Messrs.  H.  B.  Simmons,  Chief  of  the  Hydraulics  Laboratory;  F.  A. 
Herrmann,  Jr.,  Assistant  Chief  of  the  Hydraulics  Laboratory;  R.  A.  Sager 
Chief  of  the  Estuaries  Division;  and  G.  M.  Fisackerly,  Chief  of  the 
Harbor  Entrance  Branch.  Publication  of  this  report,  scheduled  for  1981, 
was  deferred  due  to  lack  of  funds.  Mr.  D.  A.  Crouse  managed  the  data 
collection  Drogram.  Dr.  R.  H.  Multer  provided  coding  for  the  digital 
filter  used.  The  computer  code  used  for  spectral  analyses  is  based  on 
a  code  written  by  Mr.  N.  W.  Scheffner.  Additional  programming,  data 
analysis,  and  preparation  of  this  report  were  performed  by  Ms.  B.  P. 
Donnell  and  Mr.  W.  H.  McAnally,  Jr. 

COL  Nelson  P.  Conover,  CE,  COL  Tilford  C.  Creel,  CE,  and  COL 
Robert  C.  Lee,  CE,  were  Commanders  and  Directors  of  WES  during  the  con¬ 
duct  of  the  study  and  the  preparation  of  this  report.  COL  Allen  F.  Grum 
USA,  was  Director  of  WES  during  the  publication  of  this  report.  Mr. 

Fred  R.  Brown  and  Dr.  Robert  W.  Whalin  were  Technical  Directors. 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 

UNITS  OF  MEASUREMENT 

Non-SI  units  of  measurement  can  be  converted  to  SI  (metric)  unit 
as  follows: 

_ Multiply _  _ By _  _ To  Obtain _ 


cubic  feet 

0.02831685 

cubic  metres 

feet 

0.3048 

metres 

SPECTRAL  ANALYSIS  OF  COLUMBIA  RIVER  ESTUARY  CURRENTS 


PART  I:  INTRODUCTION 


Objectives 


1.  The  objective  of  this  study  was  to  analyze  the  prototype  cur¬ 
rent  velocity  data  from  the  Columbia  River  estuary  and  their  relation¬ 
ships  to  local  winds.  This  segment  of  the  data  collection  effort  and 
analysis  was  to  determine  if  wind-induced  currents  within  the  estuary 
were  sufficiently  large  to  justify  their  inclusion  in  the  Corps'  efforts 
to  model  sediment  transport  in  the  estuary. 

Background 

2.  In  1976,  the  U.  S.  Army  Engineer  District,  Portland,  requested 
that  the  U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES)  conduct 
hybrid  model  study — one  using  both  physical  and  numerical  models  in  an 
integrated  manner — to  address  a  number  of  sedimentation  problems  in  the 
Columbia  River  estuary.  The  adopted  hybrid  modeling  technique  employed 
a  large  existing  physical  hydraulic  model  of  the  estuary,  a  numerical 
model  for  wave  propagation,  a  numerical  mode]  for  sediment  transport, 
and  several  analytical  techniques.  When  the  physical  processes  of  the 
estuary  were  reviewed,  it  could  not  be  established  if  wind-induced  cur¬ 
rents  played  a  significant  role  in  its  sedimentation.  It  was  decided 

to  proceed  with  the  modeling  efforts  without  attempting  to  incorporate 
wind-induced  transport,  but  to  measure  and  analyze  winds  and  currents  to 
help  define  the  relative  importance  of  winds.  Although  the  original 
list  of  study  topics  included  some  items  other  than  navigation  channel 
shoaling,  the  final  list  was  almost  exclusively  channel-related,  so  the 
goal  of  wind-current  analysis  became  limited  to  defining  the  importance 
of  wind-induced  currents  in  navigation  channel  shoaling.  This  report 
presents  results  of  the  analyses. 


PART  II:  DATA  COLLECTION  PROGRAM 


3.  Tidal  elevations  and  current  velocities  were  measured  continu 
ously  from  9  March  1978  through  6  April  1978.  Tidal  elevations  were 
measured  at  13  stations  as  shown  in  Figure  1.  Current  speed  and  direc¬ 
tion  were  measured  at  12  locations,  as  shown  in  Figure  1,  at  1,  2,  or  3 
points  in  the  water  column.  Table  1  lists  the  measurement  points,  thei 
locations,  depths,  and  starting  times.  A  plot  of  a  representative  tide 
station.  Point  Adams,  gage  3,  beginning  at  the  start  of  the  survey  is 
shown  in  Figure  2.  A  complete  presentation  of  tidal  data  is  given  by 
McAnally  and  Donnell  (in  preparation).  The  mean  freshwater  discharge 
for  the  period  of  the  survey  ranged  from  126  to  271  thousand  cubic  feet 
per  second.* 

4.  Currents  were  measured  with  Environmental  Devices  Corporation 
(ENDECO)  model  105  recording  current  meters.  The  neutrally  buoyant 
meters  were  tethered  to  a  fixed  mooring  line  by  a  5-ft-long  rope  and 
swivel  assembly.  Current  speed  and  direction  were  recorded  photograph! 
cally  at  30-min  intervals  by  an  internal  mechanism.  Current  speed  is 
sensed  by  the  rotation  of  a  horizontal  axis  propeller.  A  complete 
description  ot  the  data  collection  program  and  presentation  of  the  data 
are  given  by  McAnally  and  Donnell  (in  preparation) . 

5.  Wind  data  were  obtained  at  3-hr  intervals  for  Astoria  Airport 
Wind  speed  is  given  in  knots  and  direction  is  in  tens  of  degrees  clock¬ 
wise  from  true  north.  A  complete  presentation  of  the  data  is  given  by 
the  National  Oceanic  and  Atmospheric  Administration  (NOAA  1978). 


*  A  table  of  factors  for  converting  non-Sl  units  of  measurement  to  SI 
(metric)  units  is  presented  on  page  3. 


PART  III:  DATA  ANALYSIS  TECHNIQUES 


6.  This  portion  of  the  report  describes  cross-spectral  analysis 
and  filtering  techniques  that  have  been  applied  to  the  data  as  aids  to 
interpretation.  Spectral  analysis  is  a  complex  mathematical  technique, 
rooted  in  Fourier  analysis,  that  is  used  to  analyze  time  series  data. 

The  following  paragraphs  briefly  describe  the  analysis  technques  and 
their  interpretation.  More  detailed  treatments  can  be  found  in  Blackman 
and  Tukey  (1958),  Jenkins  and  Watts  (1969),  Bath  (1974),  and  Newland 
(1975). 

7.  This  analysis  treats  time  series  data,  a  sequence  of  data 
points  over  a  period  of  time,  as  a  general  function  that  can  be  described 
both  in  the  time  domain  and  the  frequency  domain.  This  dual  representa¬ 
tion  is  illustrated  by  the  simple  time  domain  function 

g(t)  =  C1  sin  (c^t) 

which  can  be  expressed  in  the  frequency  domain  as 

G(o)  =  for  0  = 

=  0  for  a  ^ 

This  duality  is  used  extensively  in  the  following  paragraphs. 

Fourier  Analysis 


8.  Any  function  f(t)  which  satisfies  the  following  conditions, 

can  be  represented  by  a  Fourier  series. 

a.  The  function  is  periodic  with  period  T^  such  that  the 
following  holds  for  all  time. 

f(t)  =  f(t  +  Tq) 
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Since  it  is  physically  impossible  to  acquire  a  signal 
which  begins  at  minus  infinity  and  ends  at  plus  infinity, 
it  will  suffice  to  have  periodic  data  during  an  observable 
interval  by  using  the  finite  Fourier  series. 

_b.  The  function  has  a  finite  number  of  discontinuities  in 
any  period. 

c^.  The  function  contains  a  finite  number  of  minima  and  maxima 
during  any  period. 

d.  The  function  must  be  absolutely  integrable  in  any  period. 
Otherwise  stated,  the  absolute  value  of  the  area  under 
the  curve  for  any  period  must  be  less  than  infinity. 

The  last  three  conditions  are  called  Dirichlet  conditions  (Smith  and 
Thornhill  1979).  If  f(t)  meets  these  conditions  then  it  can  be 
represented  by  a  continuous  Fourier  series  of  the  form. 


oo 


where 


n  =  0,  1 ,  2 ,  . .  . 


=  12  3 


where 

Tq  =  period 

f(t)  =  approximation  of  f(t) 
t  =  time 

9.  The  complex  form  of  the  continuous  Fourier  series  (Hamming 


1973)  is  given  by 


Equation  3a  for  is  also  known  as  the  complex  amplitude  of  the  k 


th 


component  and  is  sometimes  denoted  by  F(k) 
conjugate,  is  denoted  by  F(k)  . 


Equation  3b,  the  complex 


\  =  i\  f(t  )  sin  k  =  0,  1,...,  N-l 

p=0  0 

pT 

t  =  ~z — ~  p  =  0,  1,  2, . .  . ,  2 N-l 

p  2N  1 

2N  =  number  of  sample  points 

11.  The  continuous  Fourier  series  is  designated  by  lower  case 
coefficients  and  the  finite  form  by  upper  case  coefficients.  The 
finite  Fourier  series  coefficients  can  be  written  in  terms  of  the  con¬ 
tinuous  Fourier  series  coefficients  (Hamming  1973),  as  given  below 


A 

o 


a 

o 


+ 


CO 


m=l 


\  ~  ak  +  2  <a2Nm-k  +  a2Nm+k)  (5 

m=l 

00 

Bk  =  bk  +  2  (~b2Nm-k  +  b2Nm+k) 
m=l 

where  2N  equals  the  number  of  sample  points.  This  comparison  reveals 
that  various  frequencies  present  in  the  original  signal  are  added 
together  due  to  the  sampling.  This  phenomenon,  known  as  aliasing,  will 
be  discussed  in  paragraphs  34  and  35. 


Spectral  Analysis 

12.  The  time  and  frequency  domains  duality  is  achieved  in  Fourier 
analysis  by  use  of  the  Fourier  transform  pair.  For  the  general, 
aperiodic  function  f(t)  ,  the  general  transform  pair  is 


11 


/2tt  ^ 


F(a)  exp  (iot)  da 


■±/ 


f(t)  exp  (-iat)  dt 


where  a  equals  frequency. 

13.  F(a)  and  f(t)  may  be  thought  of  as  the  frequency  domain 
and  the  time  domain  representations  of  a  general  function  F  .  They 
represent  the  same  function  expressed  in  different  domains,  and  Equa¬ 
tions  6a  and  6b  permit  transition  from  one  domain  to  the  other. 

Power  density  spectra 

14.  The  frequency  domain  function  F(o)  in  Equation  6  is  the 
amplitude  spectrum  of  f(t)  .  For  an  f(t)  composed  of  purely  periodic 
functions,  F(a)  will  consist  of  a  series  of  lines  representing  the 
amplitude  of  the  various  frequency  components  as  shown  in  Figure  3. 


FREQUENCY 


Figure  3.  Amplitude  spectrum  of  a  sum  of  purely 
periodic  functions 
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15.  The  phase  spectrum  of  a  harmonic  f(t)  is  given  by  the  fol¬ 
lowing  equation  and  represents  the  phase  shift  of  each  of  the  periodic 
components  from  the  beginning  of  the  record 

c(ok)  =  (7> 

and  its  plot  is  similar  to  that  for  F(o)  . 

16.  The  auto-spectral  density  function  (ASD) ,  often  called  the 
"power  density,"  expressed  as 


is  a  useful  measure  of  the  distribution  of  energy  over  the  various 
component  frequencies.  Phase  information  is  not  retained.  The  term 
power  spectral  density  is  a  holdover  from  communication  engineering, 
where  these  techniques  were  developed.  In  most  realistic  applications 
the  spectrum  is  described  by  a  curve  of  power  density  or  energy,  leading 
to  equivalence  of  the  terms  power  spectral  density  and  energy  spectra. 
The  computations  described  herein  yield  a  power  spectral  density  func¬ 
tion,  which  when  multiplied  by  appropriate  physical  constants  would 
yield  the  power  density.  Since  the  absolute  value  of  auto  (power) 
density  is  not  a  desired  output,  we  will  use  the  auto  (power)  density 
function. 

Auto- co variance 

17.  The  auto-covariance  represents  the  lagged  product  of  a  time 
series  X(t)  (with  the  mean  removed)  multiplied  by  itself.  It  is  de¬ 
fined  as 


XX 


CO 

(t)  =  J  [X(t)  -  x][x(t  +  T)  -  x]  dt 


(9) 


where 


t  =  time  lag,  or  the  time  interval  which  the  series  is  shifted 


X  =  mean  of  the  X(t)  series 


18.  The  auto-covariance  is  used  to  detect  periodicities  in  a 
series,  because  its  value  will  be  large  if  the  set  of  observations  is 
multiplied  by  the  same  observations  one  period  apart. 

Cross-covariance 

19.  The  cross-covariance  is  the  lagged  product  of  a  time  series 
X(t)  multiplied  by  another  series  Y(t)  as  the  Y  series  is  shifted 
in  time  (lagged).  The  cross-variance  is  defined  as 


00 

-/&(, 


-  x][Y(t  +  x)  -  y] 


where  Y  equals  the  mean  of  the  Y(t)  series. 

20.  The  cross-covariance  will  exhibit  maxima  at  periods  that  are 
common  to  both  series. 


Auto-spectral  densit} 


21.  The  auto-spectral  density  (ASD)  of  a  function  is  the  frequency 


representation  (Fourier  transform)  of  the  auto-covariance,  defined  as 


*xx<0) 


<J>xx(t)  exP  (-i0T)  dT 


(11a) 


22.  The  auto  spectrum  does  not  retain  phase  information.  If  the 
frequencies  in  the  original  record  are  multiples  of  each  other,  the 
function  will  have  a  maximum  at  each  time  lag  that  corresponds  to  one 
or  more  periods  of  the  highest  harmonic. 

23.  The  auto  spectrum,  a  real  function,  can  also  be  expressed  as 

$xx(o)  =  Vo)  Fx(a)  =  lFx(o)i2  (llb) 

where  F  (o)  equals  the  complex  amplitude  spectrum  for  the  X  series. 

A 

Cross-spectral  density 

24.  The  cross-spectral  density  (CSD)  is  the  frequency  domain 


representation  (Fourier  transform)  of  the  cross-covariance  and  is 
defined  as 


$XY(0) 


(x) 


exp 


(-iox)  dx 


(12) 


Another  representation  of  the  CSD  is 


$XY(0)  =  FX(o)FY(o) 


(12b) 


25.  The  cross  spectrum  will  exhibit  maxima  at  frequencies  common 
to  both  series.  It  is  an  expression  of  energy  in  the  mutual  frequency 
components  of  X  and  Y  .  The  cross  spectrum  is  a  complex  function 
with  a  real  part,  the  co-spectrum,  and  an  imaginary  part,  the  quad- 
spectrum.  The  co-spectrum  represents  the  in-phase  average  product  of 

X  and  Y  .  It  will  dominate  (absolutely)  the  cross  spectrum  of  two 
functions  that  are  close  to  the  same  phase  in  time.  The  quad-spectrum 
represents  the  out-of-phase  average  product  of  X  and  Y  . 

Phase  spectra 

26.  The  phase  spectrum,  given  by 


£xy(0)  = 


tan 


■l[ 


quad-spectrum 

co-spectrum 


(o)~ 

(a)  _ 


(13) 


is  the  lag  between  components  of  the  same  frequency  in  X  and  Y  .  If 
X  and  Y  have  no  frequency  in  common,  e  (o)  will  be  constant  between 
a  =  — tt / 2  and  tt/2  .  The  response  function  amplitude  (RFA)  is  given  by 


l*XY(o) 

^(0) 


(14) 


For  a  linear  system  in  which  an  input  X(t)  produces  an  output  Y(t)  , 
knowing  X(t)  ,  the  RFA,  and  the  phase  spectrum  permits  computation  of 


Y(t)  .  The  Y(t)  spectrum  is  the  product  of  the  X(t)  and  the  RFA 
spectra.  The  phase  spectrum  represents  the  phase  shift  of  each  fre¬ 
quency  component  with  respect  to  the  input. 

Coherency 

27.  The  coherency  function  of  two  functions  is  defined  by  the 
coherence  squared, 


2  $XY(a) 

YXY(0)  “  *xx(o)*YY(o> 


which  has  values  between  0  and  1  .  If  X  and  Y  are  linearly  cor- 

2  2 

related  y  -*■  1  .  If  they  are  completely  uncorrelated  y  =  0  .  The 
Ai  AY 

presence  of  a  peak  at  the  same  frequency  in  the  auto-spectral  density 
curves  of  both  X  and  Y  is  insufficient  for  correlation.  Perfect 
linear  correlation,  a  coherence  function  equal  to  1,  is  achieved  when 
the  X  and  Y  series  are  approximated  by  the  same  Fourier  series 
(Equation  l)  except  for  a  constant  multiplier;  thus,  the  magnitude  of 
the  cosine  terms  with  respect  to  the  sine  terms  will  be  equal  in  the 
two  series. 

28.  In  summary,  spectral  analysis  transforms  time  series  data 
into  the  frequency  domain  and  computes  the  distribution  of  energy  over 
frequency  of  individual  data  records  (auto-spectral  density)  and  of  the 
product  of  two  data  records  (cross-spectral  density).  It  quantifies 
amplitude  and  phase  relationships  between  the  two  data  records  (response 
function  amplitude  and  phase  spectrum,  respectively)  and  quantifies  the 
degree  of  linear  correlation  between  the  records  (coherency  squared). 

Practical  Aspects  of  Spectral  Analysis 

29.  The  preceding  paragraphs,  concerned  with  the  theoretical  re¬ 
lationships  of  functions,  do  not  address  practical  aspects  of  data 
analysis.  An  actual  time  series  data  record  consisting  of  a  finite  num¬ 
ber  of  observations  at  discrete  intervals  between  data  points  introduces 
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several  problems  into  the  analysis.  This  section  addresses  these  prob¬ 
lems  and  how  they  are  handled. 

30.  Time  series  consisting  purely  of  constituent  periodic  func¬ 
tions  can  be  represented  as  line  spectra  in  which  all  of  the  energy  is 
located  at  precise  frequencies  and  the  power  spectrum  can  be  represented 
as  a  series  of  lines  on  a  [F  x(a)]2  versus  o  graph  (as  in  Figure  3). 
Aperiodic  and  periodic  functions  masked  by  noise  or  inadequate  record 
lengths  appear  as  continuous  spectra  in  which  the  spectrum  is  blurred 

and  occurs  as  a  curve  of  power  density  on  the  [Fx(0>]2  versus  o  graph. 
Noise 

31.  The  simplest  definition  of  noise  is  that  given  by  Godin  (1972) 
as  any  variation  in  f(t)  that  is  not  desired.  Noise  includes  measure¬ 
ment  error  and  actual  physical  phenomena  that  confuse  the  data.  Applica¬ 
tion  of  Godin's  definition  means  that  the  intent  of  the  analysis  deter¬ 
mines  what  is  noise  and  what  is  not.  If  the  objective  is  to  examine 
nontidal  phenomena,  then  the  tidal  oscillations  become  noise.  Since 
noise  obscures  those  data  record  features  of  interest,  it  is  desirable 

to  remove  it  by  filtering  or  to  minimize  its  impact. 

32.  Noise  occurs  in  the  data  and  in  the  spectrum  because  of  finite 
record  length,  the  lag  window's  characteristics  (paragraph  38),  and 
others.  This  noise  sometimes  leads  to  nonsensical  values  of  the  spectral 
estimates,  such  as  a  negative  auto-spectral  density,  coherency,  or  RFA. 
Although  theoretical  relationships  dictate  that  ASD  and  RFA  always  be 
positive,  and  that  coherency  squared  be  between  0  and  1,  actual  data 
analyses  often  result  in  some  values  falling  outside  these  bounds. 

33.  A  source  of  noise  in  this  report's  data  is  the  12-day  (or 
longer)  period  events  that  appear  in  the  0.0  to  0.003  frequency  range. 

It  is  considered  to  be  noise  in  the  sense  that  storm  events  that  are 
likely  to  affect  currents  within  the  estuary  will  have  a  duration  of 
less  than  12  days;  thus,  those  processes  of  12  day-periods  or  greater 
are  of  secondary  interest. 

Sampling  interval 

34.  The  time  between  successive  data  points  is  the  sampling  in¬ 
terval  At  which  limits  the  period  of  a  function  f(t)  that  can  be 
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accurately  represented  by  the  sampled  data  points.  The  minimum  period 
that  can  be  distinguished  is  2At  (maximum  frequency  l/(2At).  Smaller 
periods  (higher  frequencies)  are  aliased  as  longer  periods  because  the 
sampling  interval  makes  them  appear  longer. 

35.  Aliasing  is  a  part  of  our  everyday  world.  A  common  example 
of  aliasing  occurs  frequently  in  moving  pictures  when  fast -spinning 
wheels  are  phocographed  by  cameras  whose  shutter  speed  makes  the  time 
interval  of  sampling  (framing)  too  large  and  gives  the  illusion  of  the 
wheels  gradually  slowing  down,  stopping,  and  then  turning  backward.  The 
maximum  frequency  that  can  be  accurately  detected  by  sampling  interval 
At  is  called  the  Nyquist  frequency  or  folding  frequency 


It  is  called  the  folding  frequency  because  frequencies  higher  than  ff 
are  folded  onto  lower  frequencies  in  the  spectrum.  Thus,  energy  at  odd 
integer  multiples  of  ff  will  be  aliased  as  ff  and  energy  at  even 
multiples  will  be  aliased  at  zero  frequency  as  shown  in  Figure  4. 

Equation  15  says  that  we  must  have  at  least  two  samples  in  the  shortest 
period  present  in  order  to  avoid  aliasing. 

Record  length 

36.  For  a  purely  periodic  function  a  finite  record  length  equal  to 
a  multiple  of  the  fundamental  periods  can  be  used  to  completely  describe 
the  function.  Aperiodic  functions  or  basically  periodic  functions  with 
nondeterministic  information  cannot  be  completely  described  by  a  finite 
record  length  since  one  segment  of  record  will  have  somewhat  different 
spectra  than  another  segment.  Thus,  spectra  of  functions  of  the  latter 
type  (most  practical  data)  are  merely  estimates  of  the  actual  spectra. 

37.  Since  periodical  record  length  is  finite,  the  maximum  lag  T 
is  also  finite  and  limited  by 


T  <  T 
m  n 

where  T  is  the  record  length.  In  practical  terms  T  must  be  much 
n  m 

smaller  than  T  with  the  classical  limitation  being 
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Figure  4.  Spectrum  folding  due  to  aliasing 
(after  Blackman  and  Tukey  1958) 
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although  this  Limit  is  often  exceeded.  Finite  record  lengths  and  small 

maximum  lags  smooth  (smear)  the  spectral  lines. 

38.  Finite  record  lengtli  means  that  the  limits  on  thi  integral  in 

Equations  9  and  10  are  reduced  from  (-«■  to  °°)  to  (-T  ,  to  T  .  -  T 

n  /  2  n  /  2  m 

This  is  equivalent  to  multiplying  the  equation  by  D(t)  where 


WR(x)  = 


1.  T  >  T 


0 ,  x  <  T 


which  is  a  boxcar  function  as  show-n  in  Figure  5  and  is  called  the  lag 
window  since  it  is  the  window  through  which  we  view  the  covariance,  and 
thus  it  limits  our  knowledge  of  the  covariance.  When  the  cross¬ 
covariance,  1 ' (t)  ,  (the  prime  indicating  that  it  is  an  apparent  or 

A  I 

approximate  1>  (x))  is  transformed  into  the  frequency  domain  ($■'  (a)), 

A  Y  X  Y 

w ( T )  is  transformed  into  W  (o)  ,  the  spectral  window,  which  relates 

K 

the  actual  cross-spectral  density  to  its  estimated  value  by 


>XY(0)  *  V^XY^ 


W  (  ’)  is  the  frequency  domain  window  through  which  we  view  the  spectrum 
and  it  >mears  the  spectrum,  spreading  the  energy  at  each  frequency  to 
adjacent  frequencies  (leakage),  adding  negative  energy  to  some  frequen¬ 
cies,  and  adding  some  high  frequency  noise  to  the  spectrum.  Figure  5 
illustrates  the  boxcar  window  as  it  appears  in  both  the  time  and  fre¬ 


quence  domains. 


The  frequency  plot  clearly  shows  the  leakage-causing 


lobes  to  either  side  of  the  pr.mary  peak. 

.  Since  we  automatically  have  a  lag  window  by  taking  a  finite 
re  ord  length,  it  would  be  to  our  advantage  to  design  a  window  that 
would  minimi/e  the  undesirable  effects  associated  with  the  window.  If 
the  boxcar  lag  window  is  used,  the  spectrum  is  called  the  raw  spectrum. 


There  are  a  number  of  other  windows  that  accomplish  various  effects. 

The  window  chosen  here  is  the  Tukey-Hann  window  shown  in  Figure  5,  which 
is  designed  to  minimize  leakage  from  strong  peaks  of  power  density.  Ap¬ 
plication  of  the  Tukey-Hann  window  results  in  a  spectrum  that  is  smoothed. 
The  lag  window  through  which  we  view  the  cross-covariance  when  applying 
the  Tukey-Hann  equation  is: 


x  <  -T 


ws(t)  =4  +  cos  (r)]  -Tm±  T  ± 


(20) 


T  >  T 


Similarly,  when  the  cross-covariance  is  transformed  into  the  frequency 
domain,  Wg(x)  is  transformed  into  Wg(x)  ,  the  Tukey-Hann  (smooth) 
spectral  window.  The  lower  half  of  Figure  5  illustrates  the  Tukey-Hann 
window  in  both  the  time  and  frequency  domain.  Note  the  decrease  in 
amplitude  of  the  peak  and  side  lobes  of  the  Tukey-Hann  window  as  compared 
with  the  boxcar  window.  A  smoothed  spectrum  provides  benefits  by  mini¬ 
mizing  spurious  oscillations,  but  carries  the  disadvantage  that  smearing 
of  energy  at  adjacent  frequencies  increases  the  effective  bandwidth,  thus 
decreasing  ability  to  resolve  energy  peaks  that  are  not  widely  separated 
in  frequency. 

40.  One  further  effect  of  record  length  deserves  mention.  Some 

aperiodic  functions,  (for  example,  a  step  function)  when  spectrally 
analyzed,  will  contribute  to  the  power  density  at  a  frequency  correspond¬ 
ing  to  1/T  Changing  the  record  length  will  cause  the  power  density 

n 

associated  with  the  function  to  move  to  a  new  location  corresponding  to 
the  new  record  length.  This  phenomenon  does  not  eliminate  the  useful¬ 
ness  of  the  spectral  analysis  results,  but  may  confuse  interpretation 
because  the  frequency  is  not  a  property  of  the  function. 

Resolution  ve rsus  stability 

41.  Resolution  is  a  measure  of  how  well  various  components  are 
defined  or  are  concentrated  at  the  correct  frequency.  Energy  density 
estimates  are  obtained  at  Aa  =  1/(2T  )  and  adjacent  estimates  are 


overlapped  considerably  for  the  useful  windows,  but  estimates  2Ao  apart 
do  not  have  serious  overlap  so  we  can  define  the  resolution  to  be 


R  =  (21) 

m 

Thus  a  larger  lag  (large  T  )  leads  to  better  resolution. 

m 

42.  Stability  is  a  measure  of  '■he  confidence  with  which  we  can 
view  the  estimates  of  energy  density.  It  is  possible  to  have  high 
resolution  but  poor  estimates  because  of  low  stability.  Physical  ex¬ 
pressions  of  stability  are  less  obvious  than  for  resolution,  but  it 
can  be  said  that  a  longer  useful  record  length  (T  -  T  )  results  in 
greater  stability.  Thus,  it  conflicts  with  the  requirement  for  resolu¬ 
tion.  The  usual  criterion  for  stability  is  given  by  Equation  17. 

43.  Since  the  maximum  lag  has  conflicting  requirements  on  its 
length  due  to  resolution  and  stability,  choosing  the  best  T^  can  often 
be  a  problem.  This  is  partly  overcome  by  window  closing  in  which  suc¬ 
cessive  spectral  computations  are  made  for  larger  and  larger  lag  periods 
until  the  best  maximum  lag  is  found. 

44.  Stability  may  be  more  quantitatively  expressed  in  terms  of 
confidence  intervals  and  variance.  Confidence  limits  define  the  proba¬ 
bility  that  the  true  value  of  the  estimated  variable  is  within  a  given 
range  of  values.  Thus,  they  provide  a  quantitative  measure  of  the 
reliability  of  the  estimation  technique.  Computation  of  confidence 
limits  requires  use  of  the  variance  function. 

45.  The  variance  erf  the  smoothed  phase  spectrum  estimate  can  be 
expressed  as  a  function  of  the  degrees  of  freedom  (v) 


who  re 


standardized  bandwidth 
record  1 engt h 


maximum  lag  time 


The  Tukev-llann  window,  alternately  called  either  Tukev  or 


..L..  ..  1.1  ,...  V- 
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Hann,  has  a  standardized  bandwidth  of  1.333.  The  boxcar  window's 
standardized  bandwidth  is  0.5. 
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47.  The  auto-spectrum  estimates  have  a  x  (chi-square)  distribu¬ 
tion  with  v  degrees  of  freedom.  There  is  a  (1  -  a)  probability  that 
the  true  value  of  the  auto-spectral  density  lies  between 


v  ^(o) 

upper  limit  =  — -z— - r 

^  -  f) 


V 

lower  limit  =  — -■/■■■.-  — 


where 

v  =  degree  of  freedom 

a  =  probability  that  the  true  value  lies  beyond  the  specified 
limits 

$X^(c!)  =  smooth  ASD  (i.e.,  ASD  using  the  Tukey-Hann  spectral  lag 
window) 

48.  The  amplitude  and  phase  of  the  respone  function  are  dis¬ 
tributed  as  F  (Fisher  distribution)  with  2  and  (v  -  2)  degrees  of 
freedom.  The  confidence  interval  limits  of  the  amplitude  are 


(23a) 


(23b) 


upper  limit  =  A  (a) 

X  i 


(24.i ) 


(24b) 
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where 


F„  „  =  Fisher  distribution  value 
2 ,  v-2 

=  smooth  response  function  amplitude 
YXY^°^  =  smoot^  coherency  squared 

The  limits  of  the  confidence  interval  on  the  true  phase  are  given  by 
Liu  (1974)  as 

upper  limit  =  e(o) 


+  arc  sin 


J 


[v-2  F2,v-2(1  "  a)] 


1  ~  yxy(o) 

YXY(o) 


>1/2 


(25a) 


lower  limit  =  e(o) 


-  arc  sin 


J 
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[v-2  F2,v-2(1  -  a)] 


vl/2 


1  “  YXY(0) 


yxy(0) 


(25b) 


where 

t(o)  =  smooth  phase  in  degrees 
Y  (o)  =  smooth  coherency  squared 

A  i 

The  expressions  for  upper  and  lower  confidence  limits  illustrate  that 
longer  record  lengths  and  shorter  maximum  lag  lengths  lead  to  a  desirable 
narrowing  of  confidence  limits  about  the  estimated  value  of  the  param¬ 
eter.  They  also  show  the  trade-off  between  resolution  and  stability  in 
choosing  a  lag  window. 

49.  The  relationship  between  resolution  and  stability  is  nicely 
demonstrated  by  a  paradox  of  choosing  a  lag  window.  A  boxcar  window, 
with  a  standardized  bandwidth  of  0.5,  permits  better  resolution  than 
the  Tukey-Hann  window  which  has  a  standardized  bandwidth  of  1.33.  The 
Tukey-Hann  window's  larger  bandwidth  means  higher  degrees  of  freedom 
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(Equation  22)  and  thus  narrower  confidence  limits  than  the  boxcar  window 

with  the  same  T  and  T  .  Narrower  confidence  limits  indicate  more 
n  m 

reliable  answers,  which  is  so,  but  the  estimates  obtained  by  the  Tukey- 
Hann  window  are  more  reliable  estimates  for  the  average  over  a  wider 
frequency  band,  which  means  that  the  estimates  may  be  worse,  not  better. 
It  is  therefore  advisable  to  examine  both  raw  and  smoothed  spectra  to 
see  if  improved  reliability  obscures  important  features  of  the  spectrum. 

Low  Pass  Filter 


50.  A  30-point  nonrecursive  digital  low  pass  filter  was  used  to 
remove  shorter  period  oscillations  from  the  long-term  current  and  wind 
data.  The  filter  is  designed  such  that  oscillations  of  periods  longer 
than  36  hr  are  passed  while  those  shorter  than  36  hr  are  removed  from 
the  record.  A  cutoff  period  of  36  hr  was  selected  so  as  to  remove  tide 
effects  from  the  record. 

51.  An  ideal  low  pass  filter  fitting  the  above  description  would 
be  a  square  wave  in  the  frequency  domain,  with  a  value  of  1  for  fre¬ 
quencies  less  than  1/36  cph  and  zero  for  frequencies  greater  than 

1/36  cph.  The  filter  used  here  approximates  a  square  wave  by  a  Fourier 
series.  The  filter  was  developed  by  use  of  Fourier  series  coefficients 
based  on  the  interval  of  0  -  it  and  adjusted  by  the  Lanczos  correction 
factor.  The  coefficients  were  derived  from  the  following  relations: 


B 


o 


(x)dT 


where 


f(T) 


B 

n 


o 


(nx )dx 


1.0  if  0£xj<tt/6 
0.0  if  XT / 6  <  x  <  TT 


(26) 
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52.  The  Lanczos  correction  factor. 


,  is  of  the  form: 
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where 

n  =  1,2,3,...  N 

N  =  half  the  number  of  data  points  used  to  apply  the  filter  to  a 
single  data  value 

The  final  form  of  the  filter,  as  applied  to  some  data  set  at  time  t  i 

30 

PF(t)  -  80P(t)  -  i  £  Bn“n[P(t-n4t)  +  P<t-h,4t)]  <28 

n=l 


where 

P  =  point  values  of  the  original  data 
Pp  =  filtered  point  values 

Figure  6  shows  the  frequency  domain  response  curve  for  the  30  point 
filter.  As  shown  by  Equation  28,  obtaining  a  filtered  data  value  in 
the  time  domain  consists  of  summing  weighted  values  of  the  raw  data 
value,  the  30  preceding  values,  and  the  30  succeeding  values.  This 
results  in  reduction  of  the  data  record  length  by  60  points,  30  on  each 
end.  Since  the  filter  is  symmetrical  about  the  data  point,  it  does  not 
cause  any  phase  shift  in  the  filtered  data. 

53.  The  filter  is  applied  to  raw  data  in  the  time  domain  by 
Equation  28.  Its  application  is  equivalent  to  multiplying  the  filter's 
frequency  domain  representation  by  the  data  set's  frequency  domain 
represcnt.it  ion . 


PART  IV:  ANALYSIS  PROCEDURES 


54.  Preparation  of  current  velocities  for  analysis  consisted  of 
the  following  steps. 

a.  Photographic  recordings  of  current  meter  sensor  readings 
from  each  meter  were  reduced  by  Environmental  Devices 
Corporation,  manufacturer  of  the  meters.  Current  speed 
(fps)  and  direction  (clockwise  degrees  from  magnetic 
north)  toward  which  the  current  was  flowing  at  30-min 
intervals  were  written  to  a  magnetic  tape  which  was  sent 
to  WES. 

jb.  The  data  were  plotted  as  time-histories  and  checked  for 
reasonableness. 

c_.  Component  velocities  in  the  east-west  and  north-south 
directions  were  computed  and  stored.  For  the  east-west 
components,  eastward  flowing  components  were  given  a  posi¬ 
tive  sign.  Northern  directions  for  the  north-south  compo¬ 
nents  were  also  denoted  positive.  Directions  west  and 
south  were  given  a  negative  sign.  North-south  component 
analyses  were  performed  only  for  those  station  locations 
for  which  a  significant  wind  fetch  length  offered  oppor¬ 
tunity  for  wind-induced  currents  in  that  direction. 

d_.  The  low  pass  filter  described  in  paragraphs  50-53  was 

applied  to  the  component  velocities  and  the  filtered  out¬ 
put  was  stored. 

55.  In  the  analysis,  the  data  sets  resulting  from  step  (a) — 
current  speed  and  direction  received  at  WES — are  referred  to  as  original 
data.  Data  sets  consisting  of  component  velocities  are  identified  by  a 
code  that  consists  of  the  letter  C  followed  by  a  two-letter  designation 
of  either  east-west  components,  EW,  or  north-south  components,  NS,  and 
ending  with  station  and  depth  designation.  Thus  a  data  set  of  east-west 
components  of  currents  at  sta  Tl,  bottom,  is  identified  by  CEWT1B.  Data 
sets  resulting  from  filtering  in  step  (d)  are  identified  by  the  same 
code  preceded  by  an  F.  The  filtered  version  of  CEWT1B  is  therefore  iden¬ 
tified  by  FCEWT1B.  Time  in  the  data  sets  is  expressed  in  hours  from  the 
starting  time  of  each  original  data  set  as  listed  in  Table  1. 

56.  Wind  velocities  at  Astoria  Airport  for  the  period  of  interest 
were  processed  by  these  steps: 

a.  Wind  speed  (knots)  and  direction  (clockwise  degrees  from 
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true  north)  from  which  the  wind  was  blowing  were  key¬ 
punched  from  copies  of  published  climatologic  records 
(NOAA  1978). 

t).  Speeds  were  converted  from  knots  to  fps  and  directions 
were  converted  to  degrees  (from  magnetic  north)  and  to 
the  direction  toward  which  the  wind  was  blowing  to  make 
them  consistent  with  current  data. 

c_.  This  process  was  identical  to  step  (c)  for  current 
velocities . 

d_.  For  each  individual  current  station,  there  existed  a  cor¬ 
responding  wind  file  with  the  same  starting  date  and  time. 
This  was  necessary  to  prevent  a  bias  phase  lag.  The 
titling  procedure  in  this  report  does  not  distinguish  each 
individual  wind  file  needed  to  correspond  with  the  various 
starting  times  and  record  lengths  of  the  current  stations, 
which  explains  slight  variations  in  the  ASD  for  the  seven 
various  wind  files.  For  example,  the  wind  file  accompany¬ 
ing  sta  4  and  5  has  a  starting  time  of  3-9-78  at  1000  hr 
while  the  wind  file  accompanying  sta  7  has  a  starting 
time  of  3-8-78  at  1600  hr.  The  title  procedure  used 
identifies  component  wind  velocity  by  the  lable  ASTEW  and 
ASTNS  regardless  of  starting  times.  Similarly,  filtered 
versions  are  identified  by  FASTEW  and  FASTNS. 

57.  The  preparation  process  described  above  resulted  in  four  data 
sets  for  each  record  of  original  data.  Then  cross-spectral  methods  were 
applied  to  various  combinations  of  two  data  sets  at  a  time  in  order  to 
determine  the  characteristics  of  each  data  set  and  the  relationships 
among  them.  The  boxcar  and  Tukey-Hann  lag  windows  were  employed. 

58.  Current  data  were  recorded  at  30-min  intervals  while  wind 
data  were  recorded  at  3-iir  intervals.  Data  points  at  3-hr  intervals 
were  used  for  both  in  the  succeeding  analyses.  A  3-hr  time  increment  re¬ 
sults  in  a  folding  frequency  of  0.167  cph. 

59.  Usable  data  record  lengths  ranged  from  486  to  507  hr.  A 

maximum  lag  time  of  150  hr  was  chosen  after  trials  of  30,  60,  90,  120, 

150,  and  180  hr.  The  150-hr  lag  provided  spectral  estimates  at  frequency 

intervals  of  0.00333  cph  between  0  and  the  folding  frequency  of  0.167  cph. 

A  lag  as  large  as  this  violates  the  10  percent  rule  of  Equation  17  but 

was  necessary  in  order  to  adequately  resolve  the  region  of  interest 

between  0  and  0.028  cph  (1/36  hr).  Units  of  ASD  and  CSD  in  these  calcu- 
2 

lations  are  (fps)  -hr,  which  is  proportional  to  energy  density.  To 


PART  V:  RESULTS  AND  DISCUSSION 


60.  Results  of  data  preparation  and  analysis  are  displayed  in 
Tables  2-26  and  are  discussed  in  the  following  paragraphs.  The  original 
current  data  records  are  listed  in  McAnally  and  Donnell  (in  preparation). 

61.  The  following  describes  trials  with  north-south  and  east-west 
components,  both  filtered  and  unfiltered,  of  original  wind  and  current 
data,  with  several  of  the  spectral  functions  presented  in  the  preceding 
section.  All  sites  were  considered  for  east-west  analysis.  Velocity 
sta  Til  was  excluded  from  E-W  analysis  because  of  cross-channel  wind 
obstructions  near  that  station.  Stations  believed  to  have  sufficient 
wind  fetch  to  permit  significant  north-south  component  wind-induced  cur¬ 
rents  were  Tl,  TIB,  T4,  T5 ,  T5B,  T7,  Til,  and  TUB.  Sta  T3  was  not 
analyzed  due  to  large  gaps  in  the  data.  In  order  to  fully  identify  a 
particular  result,  a  string  of  modifiers  indicating  the  data  processing 
steps  and  the  spectral  parameter  is  required.  Considerable  economy  of 
prose  has  been  achieved  by  use  of  the  data  set  codes  described  in  para¬ 
graphs  55  and  56  and  abbreviations  for  spectral  parameters.  Data  set 
codes  are  identified  the  first  time  they  are  used.  Abbreviations  used 
are  ASD  for  auto-spectral  density  function,  (recall  synonymous  term 
power  spectral  density  described  in  paragraph  16)  CSD  for  cross-spectral 
density  function,  RFA  for  response  function  amplitude,  and  coherency  for 
coherency  squared.  Ninety  percent  confidence  limits  for  the  smooth  RFA 
and  ASD  were  computed  and  are  presented  in  the  tables  under  the  heading 
"Conf idence. " 


Astoria  Winds 


62.  Wind  speeds  and  directions  at  Astoria  Airport  (Figure  ])  are 
shown  in  Plate  1.  The  zero  hour  corresponds  to  3-8-78  hour  L600  PST. 
file  original  wind  data  are  listed  in  McAnally  and  Donnell  (in  prepara¬ 
tion).  The  plot  reflects  the  direction  (referencing  magnetic  north)  to 
which  the  wind  was  blowing.  The  plot  of  direction  versus  time  gives  a 
description  of  the  data  after  the  180-deg  adjustment  process  discussed 
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in  paragraph  56.  The  plot  of  speed  versus  time  shows  the  highest  sus¬ 
tained  wind  speeds  to  occur  between  hours  105  and  120,  hours  350  and  395, 
and  hours  645  and  670.  The  vector  plot  of  the  wind  data  is  helpful  in 
determining  a  general  direction  and  magnitude  of  the  wind.  For  instance, 
from  hours  350  to  395  the  wind  is  varying  between  0  to  30  fp.  at  300  to 
360  deg,  and  is  blowing  a  constant  70  deg  at  about  25  fps  at  hours  110 
to  120. 

63.  Plate  2  shows  east-west  components  of  unfiltered  (ASTEW)  and 
filtered  (FASTEW)  Astoria  wind  velocities.  The  smaller  amplitude, 
smoother  curve  is  the  filtered  data.  Note  that  the  filtering  process 
cannot  be  applied  to  the  first  90  hr  and  last  90  hr  of  the  data  record 
without  causing  a  phase  shift;  therefore  the  record  is  shortened  by  180 
hr.  The  largest  east-west  wind  speeds  occurred  between  hours  110  and  120 
of  the  record  and  were  toward  the  east.  Other  significant  eastward  wind 
periods  were  between  hours  200  and  300,  and  between  hours  510  and  550. 

A  significant  westward  trend  occurs  only  between  hours  330  and  400. 

64.  A  revealing  demonstration  of  the  filtering  process  is  seen  in 
Plate  2  at  about  hours  160  and  250.  At  hour  160  a  period  of  winds  blow¬ 
ing  westward  results  in  a  filtered  data  peak  with  a  magnitude  (3  fps) 
roughly  equai.  to  the  eastward  peak  at  hour  250,  even  though  the  unfil¬ 
tered  eastward  speeds  were  much  higher  than  the  unfiltered  westward 
speeds.  The  filtered  data  set  peaks  are  about  equal  because  the  high 
eastward  velocities  were  relatively  short-lived  and  separated  by  periods 
of  west-blowing  winds. 

*  65.  North-south  wind  components  (ASTNS)  and  the  filtered  residual 

(FASTNS)  are  shown  in  Plate  3.  Major  northerly  episodes  occurred  near 
hour  100,  between  hours  325  and  460,  and  after  about  hour  560.  The 
periods  of  moderate  southward  winds  are  between  hours  140  and  200,  hours 
230  and  305,  and  hours  470  and  550. 

66.  Table  2  illustrates  raw  and  smooth  cross-spectral  analyses  for 
ASTEW  versus  FASTEW.  Comparison  of  the  ASD  for  the  two  data  sets  shows 
the  filtering  ef fee t--KASTEW  ASD  values  at  higher  frequencies  are  pro¬ 
gressively  smaller  than  those  of  ASTEW,  and  several  energy  peaks  in  ASTEW 
beyond  0.03  cph  (not  shown  in  the  table)  are  essentially  removed  from 
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FASTEW.  In  the  low-frequency  portion  of  the  spectrum,  relative  magni¬ 
tudes  of  ASTEW  and  FASTEW  ASD  oscillate  in  the  pattern  of  the  filter's 
response  function  as  shown  in  Figure  5. 

67.  Examination  of  the  smooth  RFA  and  coherency  illustrates  fur- 
thet  the  filter's  effect.  The  RFA  have  values  of  about  1.0  from  0.0  to 
0.02  cph,  then  decline  rapidly  to  a  value  of  0.35  at  0.03  cph.  This 
is,  and  should  be,  the  shape  of  the  filter  as  shown  in  Figure  5,  with 
some  deviation  due  to  aliasing  and  inaccuracies  in  the  spectral  esti¬ 
mates.  The  coherencies  also  demonstrate  the  filter,  showing  a  linear 
relationship  between  ASTEW  and  FASTEW  at  low  frequencies  and  a  declin¬ 
ing  coherency  at  0.03  cph. 

68.  The  data  set  FASTEW  is  used  repeatedly  in  analyses  with  cur¬ 
rent  data,  and  its  ASD,  as  well  as  that  of  FASTNS,  is  repeated  in  a  num¬ 
ber  of  subsequent  tables  to  facilitate  comparison  with  those  of  the  cur¬ 
rents  (see  paragraph  54d  for  further  information).  At  this  point  it  is 
useful  to  examine  the  FASTEW  and  FASTNS  ASD  to  learn  what  we  may  expect 
in  the  following  analyses.  1 Vo  prominent  peaks  appear  in  the  FASTEW 
ASD--a  large  one  between  0.0033  cph  and  0.01  cph  and  a  smaller  one  at 
0.02  cph.  In  the  raw  spectrum,  relatively  little  energy  appears  at  0.0 
cph,  which  will  contain  contributions  from  periods  greater  than  300  hr. 

69.  Cross-spectral  results  for  ASTNS  versus  FASTNS  are  shown  in 
Table  3.  The  raw  FASTNS  ASD  differs  from  that  of  FASTEW  in  that  there 
is  only  one  peak  instead  of  two.  The  peak  at  0.0033  cph  is  somewhat 
narrower  and  two  times  greater  in  magnitude  in  FASTNS.  In  the  smooth 
ASD,  FASTNS  appears  to  be  very  similar  to  that  of  the  raw  ASD  except  for 
a  broadened  and  compressed  appearance  caused  by  the  smearing  effect  of 
the  Tukey-Hann  window. 

70.  The  purposes  of  this  analysis  are  to  determine  if  a  signifi¬ 

cant  relationship  between  wind  and  currents  existed,  and  if  so,  the 
magnitude  of  current  response  to  wind  forcing.  Evidence  of  such  a  rela¬ 
tionship  at  a  given  frequency  will  consist  of  a  coherency  near  1.0  (but  j 

less  than  1.1),  a  negative  phase  lag  from  wind  to  currents  of  signifi-  < 

cant  magnitude  to  allow  the  current  a  reasonable  response  time,  and  an  J 

RFA  that  is  consistent  with  known  characteristics  of  wind-induced 
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currents,  i.e.,  a  near  surface  speed  no  more  than  about  4  percent  of 
wind  speed. 

Station  T1 

71.  Original  c  rrent  data  for  sta  XI  (surface)  are  shown  in 
Plate  4.  The  predominant  ebb  flow  direction  for  the  period  is  toward 
the  west-northwest,  about  300  deg,  and  predominant  flood  flow  direction 
is  to  the  east-southeast  at  about  100  deg.  Maximum  current  speeds  were 
i.jout  5.5  fps. 

72.  Plates  5  and  6  illustrate  the  unfiltered  and  filtered  compo¬ 
nents  of  current  velocity  in  the  east-west  and  north-south  directions, 
respectively.  Lt  can  be  seen  that  the  great er-than-36-hr-per iod  cur¬ 
rents  are  a  small  fraction  of  the  measured  current  speeds.  The  filtered 
east-west  (FCEWT1)  currents  appear  to  be  approximately  balanced  whereas 
the  north-south  (FCNST1)  filtered  currents  are  predominantly  northward. 

73.  Cross-spectral  analysis  between  data  sets  FCEWT1  and  FCEWT10 
was  performed  in  an  attempt  to  identify  those  portions  of  the  spectrum 
at  Tl  which  were  attributable  to  river  runoff.  If  it  can  be  assumed 
that  (a)  the  spectrum  at  T10  is  much  less  affected  by  wind  than  that  at 
Tl,  and  (b)  the  river  runoff  portion  of  the  spectrum  should  contain  more 
energy  at  T10  than  at  Tl,  then  comparisons  of  the  spectra  should  reveal 
whether  the  energy  at  a  given  frequency  is  of  runoff  origin.  Assumption 
(a)  is  based  on  the  observations  that  T10  currents  should  be  relatively 
immune  to  north-south  wind  effects  because  the  river  is  so  narrow  at 
that  point;  that  winds  from  the  east  will  be  operating  over  a  much 
shorter  fetch  at  T10  than  at  Tl;  and  that  currents  induced  by  west  winds 
will  tend  to  be  dissipated  by  the  geometry  and  shallows  downstream  from 
110.  Assumption  (b)  is  valid  because  the  estuary  is  narrower  at  sta  T10 
and  the  same  freshwater  flow  will  create  higher  velocities  there  than  at 
sla  'll  unless  density  stratification  concentrates  seaward  flow  in  a  thin 
upper  layer.  The  data  do  not  suggest  that  the  latter  occurs. 

74.  Table  4  shows  results  of  cross-spectral  analysis  for  data 
sets  FCEWT1  and  FCEWT10.  Sta  T10  exhibits  raw  energy  peaks  between 


0  and  0.0033  cph  and  at  0.010  cph,  but  the  raw  and  smooth  auto  density 
magnitudes  are  consistently  smaller  than  those  at  correspond i ng  poinLs 
in  the  FCEWT1  spectrum.  In  the  smooth  spectra  the  ASD's  of  the  two 
data  sets  are  essentially  equal  at  0.0  cph.  This,  in  combination  with 
assumptions  (a)  and  tb) ,  suggests  that  energy  due  to  ri.er  runoff,  if 
present,  is  limited  :  o  :  re;:.,  n.  ies  well  below  (i.  0033  cph. 

t  o  •  ilesu  .  t  s  *  e  ,  r . 1  >  ‘I  —  s p ,  Ira.  a:,  a  ;  vs.  s  ^  >  i  :  ,\ t  1  ; . ana  !  (.<  Pl\  1  1  art.* 
shown  in  fable  >.  i  ae  energy  t  t'iK  erl  rations  ip.  ITFUTl  at 

about  the  same  f  requeti-  ies  as  ’  in  w  i  mi  data  except  that  the  raw  estimate.; 
exhibit  three  separate  Fti.dll  peaks  at  U.0,  0.01,  and  O.Od  cph  as  con¬ 
trasted  to  the  primary  and  secondary  peaks  which  occur  at  0.003  and 
0.02  cph  for  the  wind. 

76.  The  strongest  coherency  between  FASTFW  and  FCEWTL  occurs  at 
0.01,  0.02,  and  0.023  cph  where  the  raw  coherency  varies  between  0.75 
and  1.08.  The  raw  response  function  amplitude  (RFA)  at  those  frequencies 
is  about  0.05.  At  0.01  and  0.023  cph  the  phase  lag  is  p<  sitive  which 
indicates  that  the  current  preceded  the  wind.  However,  at  0.02  cph  the 
phase  is  negative  for  both  the  raw  and  smooth  spectrum,  indicating  that 
the  wind  preceded  the  current  by  2  to  5  hr.  Notably,  the  RFA  f or  the 
raw  and  smooth  case  at  0.02  cph  is  0.05  and  0.03,  respectfully.  Such 
amplitudes  coincide  with  observations  that  wind-induced  surface  currents 
have  a  speed  of  2  to  5  percent  of  wind  speed. 

77.  Table  6  shows  raw  and  smooth  cross-spectral  analyses,  respec¬ 
tively,  for  FASTNS  and  FCNST1.  The  current  exhibits  a  series  of  ASD 
peaks  in  the  raw  spectrum,  whereas  in  the  smoothed  spectrum  tile  peaks 
blend  into  a  continuous  curve.  Significant  raw  coherencies  (0.66  and 
0.67)  with  corresponding  negative  phase  lags  (-7  and  -1  hr)  occur  at 
frequencies  0.006  and  0.016  cph,  respectively.  Smooth  estimates  at 
these  frequencies  retain  their  negative  phase  but  coherencies  drop  con¬ 
siderably  (0.27  and  0.17). 


Station  TIB 


78.  Original  data  for  sta  TIB  (near  bottom)  are  shown  in  Plate  7. 


Flood  phase  currents  are  to  the  east  southeast,  about  100  deg,  as  at  the 
near-surface  sta  11.  Kbb  currents,  however,  are  more  to  the  west  than 
the  surface  currents.  Current  speeds  are  significantly  lower  at  the 
bottom  position,  with  only  a  few  readings  exceeding  3.5  fps.  A  period 
of  no  speed  (less  than  0.2  fps)  readings  occurred  between  hours  485  and 
530  and  the  data  are  seen  to  demonstrate  unusual  behavior  for  some  time 
prior  to  and  following  the  gap,  indicating  meter  fouling  during  the 
period . 

79.  Fast-west  components  of  TIB  currents  for  the  unfiltered 
(CFWT1.B)  and  filtered  (FCEWT1B)  are  shown  in  Plate  8.  Filtered  currents 
are  of  low  magnitude  and  exhibit  no  obvious  predominance  toward  the  east 
or  west. 

SO.  North-south  components  of  TIB  are  shown  in  Plate  9.  The  un¬ 
filtered  (CNST1B)  current  speeds  are  substantially  lower  than  east-west 
components  of  the  previous  plate,  but  the  filtered  results  (FCNS11B)  are 
roughly  the  same  magnitude  as  FCFWT1B  because  CNST1B  is  less  uniformly 
periodic  than  CEWT1B. 

81.  Table  7  contains  raw  and  smooth  spectra  for  FASTF.W  versus 
FCEW'l'lB.  FCEWT1B  demonstrates  three  ASD  peaks  at  0.0,  0.006,  and 
0.016  cpii.  Ihe  first  two  straddle  the  major  peak  at  0.003  cph  in  FASTEW 
The  third,  and  primary  peak,  coincides  with  the  wind's  secondary  peak. 

At  the  O.Olb-cph  peak  there  is  a  relatively,  high  raw  coherency  of  0.57 
and  a  negative  phase  of  8  hr.  However,  an  unusually  large  raw  UFA  of 
0.11  drops  to  0.06  when  smoothed.  The  smoothing  process  causes  the  nar¬ 
rowly  separated  peaks  to  merge,  broadening  and  flattening  the  peak  and 
reducing  the  coherency  to  0.26. 

82.  Spectral  results  for  filtered  north-south  components  of  wind 
and  sta  TIB  currents  are  shown  in  lable  8  for  the  raw  and  smooth  cases. 
The  raw  FOIST  LB  ASD  exhibits  peaks  at  0.067  and  0.016  cph.  The 
0.067-cph  frequency  reveals  not  on  y  a  significant  raw  coherency  of 
0.87  but  also  a  negative  phase  lag  (28  hr)  and  a  RFA  of  0.02.  The  roller 
encv  drops  considerably  when  smoothed.  The  fact  that  coherency  does  not 
have  to  accompany  an  ASD  peak  to  be  significant  is  demonstrated  at  the 
0.01-cph  frequency.  Here  the  raw  estimates  have  a  coherency  of  0.50,  a 


negative  phase  of  13  hr,  and  a  RFA  of  0.03.  The  smooth  estimates  main¬ 
tain  about  the  same  values  of  coherency,  negative  phase,  and  RFA. 

Station  T2 

83.  Original  data  for  sta  T2  currents  are  shown  in  Plate  10. 

Flood  currents  were  directed  to  the  east-southeast  and  ebb  currents  to 
the  west-northwest.  The  current  rose  demonstrates  an  apparent  blind  spot 
in  the  direction  data  at  about  290  deg,  which  probably  resulted  from  a 
malfunction  of  the  meter.  The  spring-neap  tidal  range  variation  and 
diurnal  inequality  are  quite  noticeable  in  the  plot  of  current  speeds — 
where  the  neap  phase  of  the  tides  clearly  occurs  at  about  hours  190  and 
570,  and  the  maximum  ebb  speeds  between  higher  high  water  and  lower  low 
water  are  twice  as  large  as  the  maximum  flood  speeds  and  the  ebb  between 
lower  high  water  and  higher  low  water.  Maximum  ebb  current  speeds  at 
sta  T2  were  less  than  5  fps,  while  maximum  flood  currents  were  less  than 
4  fps. 

84.  East-west  components  of  flow  are  shown  in  Plate  11.  Filtered 
east-west  components  (FCEWT2)  are  all  negative,  indicating  a  complete 
ebb  bias  in  the  nontidal  current  velocities  for  the  period.  Table  9  con¬ 
tains  the  spectral  estimates  for  FASTEW  versus  FCEWT2.  The  raw  ASD  for 
FCEWT2  shows  a  primary  peak  between  0.0  and  0.003  cph  and  a  smaller  peak 
at  0.020  cph.  Strong  raw  coherencies  at  0.0067  and  0.02  cph  are  accom¬ 
panied  by  positive  phase  lag.  Smoothing  results  in  negative  phase  and 
marginal  coherencies  of  around  0.45  for  frequencies  of  0.013  and 

0.016  cph. 

Station  _T2B 

83.  Original  data  for  sta  i'2B  (near  bottom)  currents  are  shown 
in  Plate  12.  FLood  currents  were  directed  slightly  south  of  east  at  a 
maximum  of  2.73  fps  and  ebb  currents  were  toward  the  west  to  west- 
northwest  at  a  maximum  of  2.75  fps.  The  speed  versus  time  graph  reveals 
a  zero  speed  between  hours  70  and  90. 


86.  East-west  components  of  the  current  are  shown  in  Plate  13. 


Filtered  east-west  components  (FCEWT2B)  are  a  mixture  of  flood  and  ebb; 
however,  there  are  three  rather  prolonged  ebb  periods  but  none  of  flood 

87.  Table  10  contains  the  spectral  estimates  for  FASTEW  and 
FCEWT2B.  Both  the  wind  and  current  have  their  raw  secondary  peaks  near 
0.023  cph.  However,  only  at  the  primary  peak  for  FCEWT2B  (0.01  cph)  is 
there  a  significant  coherency  with  a  negative  phase  lag.  However,  the 
smooth  coherency  drops  to  0.31  as  opposed  to  the  0.87  raw  estimate. 

Stations  T3,  T3C,  and  T3B 


88.  Records  at  all  three  depths  of  location  T3  were  marred  by 
gaps  in  the  data  to  an  extent  that  prohibited  spectral  analysis.  The 
original  data  for  T3,  T3C,  and  T3B  are  shown  in  Plates  14,  15,  and  16. 

Station  T4 


89.  The  original  measured  data  for  sta  T4  are  shown  in  Plate  17. 
The  ebb  flow  direction  for  the  period  is  toward  the  west-southwest  at 
about  250  deg  while  the  flood  flow  is  to  the  northeast  at  about  45  deg. 
Maximum  ebb  current  speeds  are  about  4.7  fps,  and  maximum  flood  currents 
are  about  3.2  fps. 

90.  Plates  18  and  19  illustrate  the  unfiltered  and  filtered  compo¬ 
nents  of  current  velocity  in  the  east-west  and  north-south  directions, 
respectively.  The  filtered  east-west  (FCEWT4)  currents  are  all  negative 
which  indicates  a  westward  (ebb)  flow  while  the  filtered  north-south 
(FCNt>T4)  currents  are  essentially  balanced. 

91.  Results  of  spectral  analysis  of  FASTEW  and  FCF.WT4  are  shown 
in  Table  11.  The  raw  ASH  show  energy  concentrations  in  the  two  data 
sets  Lo  be  at  frequencies  of  0.003,  0.006,  and  0.02  cph.  Estimates  with 
significant  coherency  values  have  positive  phase  values  which  implies  tin 
unlikely  event  that  the  current  preceded  the  wind,  and  are  omitted  from 

f tirther  cons  iderat  ion  . 

92.  Results  of  spectral  analysis  of  FASTNS  and  FCNST4  are  shown 


in  Table  12.  The  wind  and  current's  raw  ASD  indicate  a  large  concentra¬ 
tion  of  energy  at  the  frequency  0.003  cph.  A  large  raw  coherency  of 
0.99  is  slightly  reduced  to  0.73  when  smoothed,  but  a  positive  phase  lag 
eliminates  it  from  consideration.  Furthermore,  this  frequency  is  not  of 
interest,  as  discussed  in  paragraph  33.  Interestingly,  at  the  G.013 
frequency  there  is  a  0.48  raw  coherency  with  a  negative  lag  of  9  hr  and 
an  RFA  of  0.03,  but  the  smoothing  process  drops  the  coherency  and  RFA  to 
0.0  and  raises  the  phase  lag  to  negative  17  hr. 

Station  T5 

93.  Original  current  data  for  sta  T5,  surface,  are  shown  in 
Plate  20.  Flood  currents  were  toward  the  northeast  with  maximum  speeds 
of  about  2  fps,  and  ebb  flows  were  to  the  west-southwest  with  maxima  of 
about  3  fps.  The  smaller  neap  phase  velocities  beginning  at  about  150  hr 
noted  at  sta  T2  appear  in  the  sta  T5  record  also;  however,  sta  T5  direc¬ 
tions  for  that  period  are  shown  to  be  consistently  in  the  flood  direction 
and  the  speeds  do  not  exhibit  the  large  maximum  ebb  values  observed  at 
sta  T2.  A  long  period  of  consistent  flood  velocities  at  this  station  is 
unlikely,  so  the  data  suggest  rather  a  stuck  direction  indicator  or  foul¬ 
ing  of  the  meter,  perhaps  by  debris  or  the  mooring  lines.  Plate  21  illus¬ 
trates  the  east-west  components  of  T5  currents,  both  unfiltered  and  fil¬ 
tered  (FCEWT5).  Other  than  the  suspicious  period  of  all  flood  flows,  the 
predominant  direction  of  FCEWT5  is  ebb  (negative).  The  combination  of 
constant  flood  flows  in  one  portion  and  predominant  ebb  flows  in  the 
remainder  imposes  an  apparent  large  amplitude  oscillation  with  a  period 
slightly  shorter  than  the  length  of  the  filtered  record.  This  will  show 
up  as  a  high  energy  peak  in  ASD  at  about  0.002  cph. 

94.  North-south  components  of  T5  currents  are  shown  in  Plate  22. 

The  prolonged  flood  flow  periods  bias  the  filtered  north-south  components 
(FCNST5)  in  the  same  way  they  affected  those  in  the  east-west  directions. 

95.  Raw  and  smooth  spectra  for  FASTEW  vs  FCEWT5  are  shown  in 
Table  13.  The  expected  large  ASD  peak  does  occur  between  0.0  and  0.0033 
cph  with  a  value  far  in  excess  of  any  previously  seen  in  the  filtered 


currents.  Raw  coherencies,  with  an  accompanying  negative  phase,  do  not 
show  significant  correlation  at  any  frequency.  Smoothed  estimates  all 
have  positive  phase  and  relatively  low  coherency. 

96.  Spectral  estimates  for  FASTNS  versus  FCNST5  are  shown  in 
Table  14.  Raw  estimates  reveal  a  strange  set  of  coherencies  and  a  phase 
lag  that  oscillates  from  positive  to  negative.  Smooth  estimates  are 
more  settled.  Smooth  estimates  that  occur  at  frequencies  0.013,  0.023, 
and  0.026  all  have  significant  smooth  coherencies  and  negative  phase  lag, 
although  there  are  no  ASD  peaks  associated  with  these  frequencies. 


Station  T5B 

97.  The  original  data  for  sta  T5B  (near  bottom)  is  shown  in 
Plate  23.  The  flood  flow  is  to  the  northeast  at  a  maximum  of  2.5  fps. 

The  ebb  flow  is  toward  the  south-southeast  about  170  deg  with  a  maximum 
speed  of  4.75  fps.  It  should  be  noted  that  this  direction  is  cross 
channel  rather  than  downstream.  The  flood  direction  is  better  aligned 
with  the  channel  (northeast).  Maximum  flood  current  was  about  2.4  fps. 

98.  The  measured  and  filtered  T5B  currents  for  the  east-west  and 
north-south  components  are  found  in  Plates  24  and  25,  respectively.  The 
filtered  east-west  current  components  (FCEWT5B)  have  a  predominant  flood 
of  less  than  0.7  fps  while  the  filtered  north-south  current  components 
(FCNST5B)  indicate  southerly  ebb  flow  of  less  than  1.5  fps. 

99.  Table  15  contains  raw  and  smooth  cross-spectral  estimates  for 
FASTEW  and  FCEWT5B.  Neither  raw  nor  smooth  estimates  satisfy  all  of  the 
coherency,  RFA,  and  phase  requirements  to  qualify  for  additional 
consideration. 

100.  The  raw  and  smooth  spectral  estimates  for  FASTNS  and  FCNST5B 
are  shown  in  Table  16.  Raw  coherency  values  are  spurious  and  all  values 
with  a  negative  phase  are  much  greater  than  one.  The  smooth  estimates  at 
frequency  0.016  have  a  coherency  of  0.52,  a  RFA  of  0.03,  and  a  negative 
phase  of  3  hr,  which  can  be  considered  a  marginal  wind/current 
correlation. 
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Station  T6 


101.  Plate  26  illustrates  original  current  data  for  sta  T6;  cur¬ 
rents  there  were  essentially  east-west  oriented,  with  maximum  flood 
speeds  of  about  2  fps  and  maximum  ebb  speeds  of  about  2.4  fps.  The  neap¬ 
spring  cycle  of  tide  ranges  is  evident  in  the  current  speed  versus  time 
plot,  but  higher  ebb  velocities  corresponding  to  the  diurnal  maximum 
range  are  not  a  striking  feature  of  the  record.  Currents  at  sta  T6 
appear  to  have  been  well  channelized  as  evidenced  by  sharp  ebb  and  flood 
reversals,  with  little  rotation  of  the  current  vector  between  phases. 

102.  East-west  components  of  sta  T6  currents  are  plotted  in 
Plate  27.  The  filtered  data  set  (FCEWT6)  is  completely  in  the  ebb 
direction  with  magnitudes  of  less  than  0.4  fps,  which  is  about  the  margin 
by  which  maximum  ebb  speeds  exceeded  maximum  flood  speeds.  Relatively 
little  periodicity  is  noticeable  in  the  filtered  data  set.  Because  the 
currents  were  so  consistently  east-west,  it  was  not  necessary  to  perform 
a  north-south  analysis. 

103.  Spectral  analysis  results  of  FASTEW  versus  FCEWT6  are  shown 
in  Table  17.  Frequencies  of  0.003  and  0.013  cph  have  significant  coher¬ 
ency  and  phase  of  0.58,  -48  hr,  and  0.82,  -9  hr,  respectfully.  Smooth 
estimates  for  these  frequencies  retain  the  negative  phase  but  coherency 
drops  to  near  zero.  Furthermore,  the  0.003-cph  frequency  is  not  of 
interest  (see  paragraph  33). 


Station  T7 


104.  The  original  current  data  (surface)  for  sta  T7  are  shown  in 
Plate  28.  The  predominant  ebb  flow  is  directed  toward  the  northwest 
(about  315  deg)  at  a  maximum  speed  of  2.1  fps  while  the  flood  flow  is 
toward  the  east-southeast  (about  110  deg)  at  a  maximum  speed  of  1.5  fps. 

105.  Plates  29  and  30  illustrate  the  measured  and  filtered  compo¬ 
nents  of  currents  in  the  east-west  and  north-south  directions,  respec¬ 
tively.  The  plots  exhibit  the  small  contribution  (less  than  0.5  fps) 
made  by  currents  with  periods  greater  than  36  hr.  The  filtered  east-west 


(FCEWT7)  current  indicates  a  small  ebb  predominance  and  the  filtered 
north-south  (FCNST7)  current  indicates  a  flood-predominant  current  of 
less  than  0.5  fps. 

106.  Table  18  shows  spectral  analysis  results  between  FAS'LhW  and 
FCEWT7.  The  estimates  at  frequency  0.02  cph  reveal  a  significant  raw 
coherency  of  0.95  with  a  corresponding  phase  lag  of  negative  3  hr.  The 
coherency  maintains  a  significant  level  of  0.50  and  a  negative  3-lir  phase 
when  smoothed.  There  is  an  ASD  peak  at  0.02  cph  for  FASTKW  hut  not  for 
FCEWT7.  Although  0.01  cph  does  not  have  an  energy  peak,  the  raw  coher¬ 
ency  is  0.57.  The  phase  lag  is  -22  hr,  and  the  RFA  is  0.01.  The 
smooched  coherency  is  not  significant,  however,  and  the  phase  lag  is 
positive. 

107.  Due  to  the  geographical  location  of  sta  T7,  analysis  was  per¬ 
formed  for  north-south  currents  versus  both  east-west  and  north-south 
winds.  The  spectral  analysis  of  FASTEW  and  FCNST7  are  shown  in  Table  19. 

A  secondary  CSD  peak  occurs  at  0.02  cph  where  a  raw  coherency  of  0.49 
and  a  phase  of  -5  hr  occur.  Smoothing  has  relatively  little  effect  on 
the  estimates  at  0.02  cph. 

108.  The  spectral  analysis  of  FASTNS  and  FCNST7  are  shown  in 
Table  20.  The  estimates  at  0.01  cph  for  both  the  raw  and  smooth  cases 
reveal  favorably  consistent  phase,  RFA,  and  coherency  squared.  The  phase 
is  approximately  -13.5  hr  and  the  coherency  ranges  between  0.58  (smoothed) 
and  0.84  (raw).  At  0.016  cph,  the  raw  phase  is  -1.0  hr  with  a  0.56  co¬ 
herency  buL  the  smoothed  estimates  show  a  positive  2-hr  phase  and  a 
minuLe  coherency  of  0.04.  There  were  no  ASD  peaks  at  frequencies  of 

ei  tiler  0.01  or  0.016  cph. 

Station  T8 

1 09.  iiie  measured  current  data  for  sta  T8  are  shown  in  Plate  31. 
i  lie  predominant  ebb  flow  direction  for  the  period  is  seen  to  be  toward 
tlie  west -southwest  at  about  260  deg,  at  a  maximum  speed  of  4.0  fps.  The 
flood  currents  are  less  than  1  fps  toward  the  east-northeast. 

110.  Plate  32  illustrates  the  measured  and  filtered  east-west 


velocities  oi  sta  T8. 


The  plot  indicates  a  decisive  ebb  direction  cur¬ 


rent  averaging  about  1.3  fps  in  the  filtered  current  (FCEWT8) .  The  domi- 
nance  of  ebb  currents  over  flood  currents  at  sta  T8  through  112  is  indic¬ 
ative  of  the  predominance  of  freshwater  discharge  over  the  tide  in  the 
restricted  upstream  portion  of  the  estuarv.  Al  any  particular  station  in 
this  area,  the  duration  and  magnitude  of  flood  current  s  decrease  as  the 
freshwater  discharge  increases.  Above  some  critical  discharge,  there  is 
no  longer  a  reversal  of  flow  direction  with  the  tide,  although  the  magni¬ 
tude  of  the  ebb  currents  will  increase  and  decrease  as  a  inaction  of  the 
tide.  This  condition  also  was  noted  downstream  aL  sta  lb  and  17,  al¬ 
though  to  a  much  lesser  degree. 

111.  The  spectra]  analysis  results  for  data  sets  FAS TEW  and  FCEWT8 
are  shown  in  Table  21.  The  raw  results  show  that  the  primary  frequence 
in  common  with  both  data  sets  is  0.003  eph,  which  as  discussed  in  para¬ 
graph  33  is  related  to  term  events  and  thus  is  not  of  interest.  The 
multiple  peaks  in  the  AS1)  of  FC11WT8  are  transformed  into  one  peak  with  a 
siightly  higher  value  when  smoothed.  A  significant  raw  coherency  of 
0.94  with  an  accompanying  negative  7-hr  lag  at  0.01  eph  is  transformed 
to  a  coherency  of  1.48  and  a  positive  3-hr  phase  lag.  Although  0.01  eph 
was  one  of  the  minor  raw  ASD  peaks  for  the  current,  it  was  not  a  peak 
for  the  smoothed  ASD. 


Station  T9 


112.  The  measured  speed  and  direction  taken  at  sta  T9  (surface) 
can  be  seen  in  Plate  33.  The  predominant  flood  flow  is  toward  the  north¬ 
east  at  45  deg  at  a  maximum  rate  of  1  fps.  The  ebb  flow  is  to  the  west 
at  about  260  deg  with  a  maximum  speed  of  2.75  fps. 

113.  As  expected,  Plate  34  describing  the  measured  (CP.WT9)  and 
filtered  (FCEWT9)  current  still  reveals  a  strong  ebb  (about  0.8  fps  aver¬ 
age)  after  the  filtering  technique  was  applied.  The  spectral  analysis, 
as  shown  in  Table  22,  does  not  indicate  a  need  to  continue  investigation. 
The  only  frequency  with  significant  energy  in  both  data  sets  is  0.0033 
which  has  a  large  positive  phase  lag  and  a  raw  coherency  of  0.89.  This 


frequency,  as  described  in  paragraph  33,  is  not  of  interest.  At  the 
0.02-cph  frequency,  the  smoothed  spectral  estimate  shows  marginally  sig- 
nigicant  coherency  of  0.49  and  a  phase  lag  of  -6  hr. 

Station  TIP 

114.  The  plots  of  the  original  data  (surface)  are  shown  in 
Plate  35.  The  polar  plot  indicates  a  flood  flow  of  about  45  deg  toward 
the  northeast  at  a  maximum  of  1.5  fps  and  ebb  currents  toward  the  south¬ 
east  (240  deg)  at  a  maximum  speed  of  3  fps.  The  measured  (FCEWT10)  and 
filtered  (FCEWT10)  currents  are  shown  in  Plate  36.  As  in  the  other  up¬ 
stream  stations,  the  filtered  data  set  shows  a  prominent  ebb  (about 

0.9  fps  average) . 

115.  Table  23  shows  results  of  spectral  analysis  of  data  sets 
FASTEW  and  FCEWT10.  The  high  CSD  values  draw  attention  to  frequencies 
between  0.0  and  0.0033,  but  these  frequencies  describe  long-term  events 
(see  paragraph  33).  A  raw  coherency  of  1.0  with  a  nositive  phase  of 

10  hr  is  at  the  0.016  frequency.  The  phase  becomes  a  negative  3  hr  with 
a  0.06  coherency  when  smoothed.  These  and  other  estimates  at  this  sta¬ 
tion  can  be  eliminated  from  consideration. 

Station  Til 

116.  Plate  37  shows  the  plots  of  the  original  current  data  for 

sta  ill  (surface).  The  flood  flow  is  less  than  1.25  fps  toward  the  south- 
southeast.  The  ebb  flow  is  almost  4  fps  toward  the  north.  In  Plate  38, 
tlie  north-south  components  of  the  measured  (CN'STll)  and  filtered  (FCNSTll) 
currents  are  shown.  The  plot  reveals  a  strong  tendency  toward  the  north, 
i.e.  ebb  flow  (average  about  1.6  fps). 

117.  The  raw  and  smooth  spectral  analysis  of  FASTNS  and  FCNSTll 
are  shown  in  Table  24.  The  raw  and  smooth  coherency  estimates  are  either 
much  greater  than  1.0  or  insignificant.  Consequently,  sta  Til  is  omitted 
from  further  consideration. 
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Station  TUB 


118.  The  near-bottom  measured  current  data  for  sta  TUB  are  given 
in  Plate  39.  The  flood  flow  is  toward  the  south  at  a  maximum  velocity 
of  1  fps,  while  the  ebb  flow  is  to  the  north-northwest  with  a  maximum  of 
3  fps. 

119.  The  measured  and  filtered  north-south  components  of  the 
current  are  shown  in  Plate  40.  The  filtered  version  (FCNST11B)  demon¬ 
strates  a  residual  ebb  flow  of  about  1.3  fps. 

120.  Table  25  contains  the  spectral  estimates  for  FASTNS  and 
FCNST11B.  Other  than  the  frequencies  describing  long-term  events,  only 
one  frequency,  0.013  cph,  has  a  negative  phase  lag  (10  hr)  with  an  accom¬ 
panying  coherency  of  0.41,  which  is  marginal.  The  smooth  estimates  for 
0.013  cph  has  a  positive  12-hr  phase  and  a  minute  coherency  of  0.01. 


Station  T12 


121.  The  original  current  data  taken  at  sta  T12  (surface)  are 
shown  in  Plate  41.  The  flood  flow  is  to  the  southeast  with  a  1.5  fps 
maximum  while  the  ebb  flow  is  to  the  northwest  at  a  maximum  of  3  fps. 

The  filtered  east-west  component  of  the  current  (FCEWT12)  as  seen  in 
Plate  42  exhibits  a  predominant  ebb  flow  (about  0.4  fps  average). 

122.  Raw  and  smooth  spectral  estimates  for  FASTEW  and  FCEWT12  are 
given  in  Table  26.  Other  than  long-term  events,  significant  coherency 
values  are  not  accompanied  by  negative  phase  lags.  This  factor  elimi¬ 
nates  this  station  from  further  investigation. 


Summary  of  Results 


123.  Table  27  lists  those  data  sets  that  exhibited  either  a  raw 
or  smooth  coherence  of  about  0.5  to  1.0,  negative  phase  lag,  and  RFA 
less  than  or  equal  to  6  percent.  Nine  wind-current  combinations  satisfy 
these  criteria  at  one  or  two  frequencies  in  either  raw  or  smooth  spectra, 
but  only  three— FASTNS-FCNST1B  at  0.01  cph,  FASTEW-FCEWT7  at  0.02  cph. 


and  FASTXS-FCNST7  at  0.01  cph — satisfy  them  in  both  raw  and  smooth 
spectra.  Seven  of  the  nine  data  sets  satisfy  the  criteria  at  0.01, 

0.0167  cph,  or  both. 

124.  Also  shown  in  Table  27  are  the  wind  and  current  speed  ampli¬ 
tudes,  j  F ( o ) j  ,  associated  with  the  raw  ASD  of  each  data  set.  These 

values  were  calculated  as  described  in  paragraph  59.  Neither  group  of 
speeds  contains  large  values,  with  the  maximum  wind  speed  of  5.6  fps 
and  maximum  current  speed  0.5  fps;  however,  a  sustained  wind  of  about 

5  fps  blowing  in  the  same  direction  for  75  hr  (frequency  0.0067  cph  in 
FASTNS)  can  have  a  significant  effect  on  water  velocity.  On  the  other 
hand,  a  wiV -induced  current  of  less  than  0.5  fps  is  unlikely  to  make  a 
major  contribution  to  instantaneous  sediment  transport  in  the  Columbia 
River  estuary,  but  as  a  residual  superposed  current  lasting  multiple 
tidal  cycles  it  could  affect  short-term  sediment  transport  patterns. 

125.  Further  analysis  might  confirm  or  cast  doubt  on  the  appear¬ 
ance  of  a  causative  relationship,  and  would  permit  an  extrapolation  to 
wind  events  other  than  those  occurring  during  the  measurement  period. 
However,  present  capabilities  do  not  offer  the  opportunity  to  carefully 
evaluate  the  contribution  of  wind-induced  currents  to  sediment  transport, 
except  at  the  large  expense  of  effort  and  funds.  At  this  point  these 
results  do  not  justify  that  effort. 

126.  From  the  analysis,  it  appears  probable  that  there  was  a 
Linear  causative  relationship  between  wind  and  currents  at  sta  TIB  and  T7 
for  frequencies  of  0.01  to  0.02  cph  during  the  period  of  measurement. 
There  is  a  significant  possibility  that  a  similar  relationship  exists  at 
the  other  stations  listed  in  Table  27. 


PART  VI:  CONCLUSIONS 


L27.  Completion  of  analysis  of  the  wind  and  currents  of  the  Colum¬ 
bia  River  estuary  leads  to  the  following  conclusions: 

a.  The  strongest  supportive  evidence  of  wind-induced  currents 
occurs  at  sta  T-l  (surface)  for  the  north-south  component 
at  a  period  of  100  hr,  and  sta  T-7  (surface)  for  the  east- 
west  component  at  200  hr,  and  the  north-south  component 
at  100  hr. 

Id.  Wind-induced  currents  generally  will  be  considerably 

less  than  0.5  fps.  Such  velocities  are  unlikely  to  make 
a  major  contribution  to  instantaneous  sediment  transport, 
but  as  a  residual  superposed  current  lasting  multiple 
tidal  cycles  short-term  sediment  transport  patterns  could 
be  affected. 

c.  Based  on  the  data  analyzed,  there  is  not  sufficient  evi¬ 
dence  to  include  wind-induced  currents  in  the  Corp's 
model  of  sediment  transport. 
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Table  1 

Station  Locations  and  Conditions 


Station 

Location 

Oregon  State  Grid 

Water 

Depth 

Meter 

Depth 

Data  Begins 

No. 

X 

Y 

ft  mllw 

ft  mllw 

PST 

T1 

1,117,318 

960,690 

52 

5 

1300 

9 

Mar 

78 

TIB 

1,117,318 

960,690 

52 

47 

1300 

9 

Mar 

78 

T2 

1,134,515 

940,811 

44 

5 

1300 

9 

Mar 

78 

T2B 

1,134,515 

940,811 

44 

39 

1300 

9 

Mar 

78 

T3 

1,140,499 
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APPENDIX  A:  NOTATION 


CSD 
I'f 
f(L) 
f  (  t  ) 

f  ( o 

F  ( k> 
F(k) 
F(  ') 

A 


k  ( t ) 
C.('3) 


n 

2N 

P(t) 

P,,(t) 

R 

RFA 


Finite  Fourier  series  coefficients  where  k  is  an  integer 
Continuous  Fourier  series  coefficients  where  n  is  an  integer 
Response  function  amplitude  of  series  X  and  Y 
Smooth  A  (a) 

Auto-spectral  density 

Finite  Fourier  series  coefficients  where  k  is  an  integer 
Continuous  Fourier  series  coefficients  where  n  is  an  integer 
Standardized  bandwidth 

Coefficient  of  complex  form  of  continuous  Fourier  series  where 
k  is  an  integer 

Cons  tant 

Cross-spectral  density 

Folding  frequency  or  Nyquist  frequency 

Function  with  respect  to  time 

The  Fourier  series  approximation  of  a  function  f(t) 

Step  function  used  for  filter 

t  h 

Complex  amplitude  of  k  component 
Complex  conjugate  of  F(k) 

Fourier  transform  of  a  function  f(t) 

Amplitude  spectrum 

Complex  amplitude  spectrum  (complex  conjugate  of  F  (a)) 

X 

Fisher  (F)  distribution  with  a  numerator  of  2  degrees  of 
freedom  and  a  denominator  with  (v-2)  degrees  of  freedom 

Time  representation  of  a  function  G 

Frequency  representation  of  a  function  G 

/-I 

In  teger 

Positive  integer 
Number  of  sample  points 
The  originaL  value  at  time  t 
The  filtered  value  at  time  t 
Resolution  in  the  frequency  domain 
Response  function  amplitude 


A1 


t 

t 

p 

T 

M 

T 

n 

T 

o 

V0) 

VT) 

w  (o) 
ws(0 

X 

Y 

a 


a 

M 

n 


'XY 

-2 

YXY 


(a) 

(a) 


Aa 

~(a) 

“(  J) 


7T 


xx([) 

XY(t> 

xx(  J 

XX ( 


XY 


.(•J) 


Time 

Finite  set  of  equally  spaced  points  in  the  finite  form  of 
Fourier  series 

Maximum  time  lag 

Record  length 

Period 

Spectral  window  in  frequency  domain  of  W  (t) 

Boxcar  lag  window  in  time  domain 

Spectral  window  in  frequency  domain  of  W  (t) 

Tukey-Hann  lag  window  in  time  domain 
Mean  of  time  series  X 
Mean  of  time  series  Y 

Probability  that  the  true  value  lies  beyond  the  specified 
limits 

Coefficients  for  low  pass  filter  where  n  is  an  integer 
Coherency  square  using  Tukey-Hann  spectral  window  (i.e.  smooth) 

Coherency  squared  of  series  X  and  Y 
Sampling  interval 
Frequency  interval 

Phase  lag  using  Boxcar  spectral  window 

Phase  lag  using  Tukey-Hann  spectral  window  (i.e.  smooth) 

Phase  spectrum 
Degrees  of  freedom 
3.14159 
F  requency 

tli 

Frequency  corresponding  to  the  k  component 

Time  lag,  or  the  time  interval  which  a  series  is  shitted 

Au  to-covar  ianc.e  of  series  X(L} 

Cross-covariance  between  series  X(L)  and  Y(t) 

Auto-spectral  density  of  series  X(t) 

Smooth  t-  (u) 

Cross-spectral  density  of  series  X ( t )  and  Y(L) 

Chi-square  distribution 

I.anczos  correction  factor  for  filter 


A2 


FILMED 


